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Abstract: We study theoretically the generation of photon pairs by 
spontaneous four-wave mixing (SFWM) in photonic crystal optical fiber. 
We show that it is possible to engineer two-photon states with specific 
spectral correlation ("entanglement") properties suitable for quantum 
information processing applications. We focus on the case exhibiting no 
spectral correlations in the two-photon component of the state, which we 
call factorability, and which allows heralding of single-photon pure-state 
wave packets without the need for spectral post filtering. We show that 
spontaneous four wave mixing exhibits a remarkable flexibility, permitting 
a wider class of two-photon states, including ultra-broadband, highly- 
anticorrelated states. 
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1. Introduction 

Quantum optical technologies require photon states with specific spectral properties. For ex- 
ample, quantum information processing using linear optics is based upon the availability of 
pure-state single-photon wavepackets (TJ. Single-photon wavepackets can be prepared using 
pair generation by means of spontaneous parametric downconversion (PDC), or by sponta- 
neous four wave mixing (SFWM) [2]. In both cases individual photons are heralded by the 
detection of their siblings. The prepared photons will not normally be in pure states unless spe- 
cial care is taken to remove all correlations in every degree of freedom of the photon pairs, i.e. 
to make the two-photon state factorable! 3 1- If this is not done, the heralded photons will be in 
mixed states, and therefore unsuitable for use in quantum logic gates, which rely on Hong-Ou- 
Mandel (HOM) interference between independent photons. Typically, photon pairs generated 
in spontaneous processes exhibit significant spectral and spatial correlations due to the energy 
and momentum conservation constraints that are typical in parametric nonlinear optics. Spatial 
correlations may be minimized by use of guided-wave configurations, such as those exploiting 
nonlinear waveguides and optical fibers I4ll5ll6l 171 l8ll9l [T0l . 

Spectral correlations are, however, more difficult to eliminate. It is possible to eliminate all 
correlations in this degree of freedom for photon pairs generated by means of PDC, as first 
shown in Ref. ifTTIl . and later extended in [3 ]. This is achieved by the method of group-velocity 
matching using a broadband pump pulse. In this paper we generalize this method to the case 
of SFWM, by engineering the group velocities of the sideband photons using photonic crystal 
fibers (PCFs). PCFs are comprised of a solid silica core surrounded by a silica cladding contain- 
ing a regular array of air holes. This leads to an exceptionally high core-cladding index contrast, 
creating strong waveguide dispersion that can be tailored for a wide-range of applications lfT2ll . 

Spontaneous four wave mixing occurs in single-mode fibers with a third-order optical non- 
linearity. In this process, two pump photons are scattered from one or two distinct pump fields 
into a pair of fields, labeled signal and idler, which are spectrally and/or polarization distinct. 
The single-mode waveguide geometry leads to the suppression of correlations in the transverse 
momentum degree of freedom, as well as suppression of mixed spectral-transverse momentum 
correlations. Thus, by imposing appropriate group velocity matching constraints on the SFWM 
process in a single-mode fiber, it becomes possible to eliminate correlations in all degrees of 
freedom, resulting in factorable two-photon states. 

Previous experimental work has explored this type of source in a number of geometries 



including standard single-mode fiber and PCF. Photon number correlations between signal and 
idler fields were first observed by pumping in the anomalous dispersion region [ 13 1, although 
there was a significant background of noise photons generated by Raman scattering of the pump 
light [14|. Spontaneous Raman scattering generates photons shifted to the red of the pump 
wavelength by up to 50 THz, corresponding to the frequency of an optical phonon in glass. 
Thermal population of the phonons gives rise to scattering over a broad frequency range, and the 
process does not require phase matching, so it occurs at all pump wavelengths. Photonic crystal 
fibers allow the signal and idler fields to be widely separated in frequency from the pump, so 
that Raman noise at the lower-frequency photon wavelength is greatly reduced EHH). Phase- 
matching of SFWM with widely separated sidebands is enabled by pumping in the normal 
dispersion region, which requires the favorable dispersion properties of strongly guided waves 

nansina. 

These experiments did not, however, address the requirement of purity for generating her- 
alded single-photon states. In particular, the photons generated in those experiments would re- 
quire tight spectral filtering to eliminate correlations between signal an idler frequencies prior 
to heralding. Using such tight filtering ensures high-visibility HOM interference between pho- 
tons from separate sources, but greatly reduces the count rates. To avoid the need for filter- 
ing, the SFWM process must produce a state whose two-photon component can be written as 
the product of the signal state and the idler state, i.e., | y/") = |0) J |0) I - + ff|l) 4 |l),-, in the pho- 
ton number basis with s, i indicating signal and idler modes respectively. We call this property 
"factorability," which corresponds to the absence of correlations between the frequencies (and 
momenta) of the idler and signal photons. In physical terms, factorability implies that no infor- 
mation about, say, the idler photon (apart from its existence) can be extracted from the detec- 
tion of the signal photon, or vice-versa. To date, such factorable states have been produced only 
in PDC at a single wavelength using a particular crystal having special dispersion properties 
|[l"8l . and other techniques could in principle be used to extend possible operation wavelengths 
J4] [19] |20l HH [22 23 1 . We show here that factorability can be achieved at a much wider range 
of wavelengths by using PCF. 

In addition to factorability, PCF enables the production of a wide range of spectral cor- 
relations in the two-photon component of the state. Two extremes are possible, a spectrally 
correlated state and a spectrally anti-correlated state. The former is a resource for quantum- 
enhanced quantum positioning [24|, whereas the latter is of importance for applications relying 
on time-of-arrival differences between two optical modes, such as optical coherence tomogra- 
phy ESI . We will show that the anti-correlated case can be made ultra-broadband by tailoring 
the higher-order dispersion. Thus we see that the ability to engineer the dispersion in PCF leads 
to a flexible system that can generate states in two widely disparate regimes (factorable and 
highly correlated), and indeed in a very general class of intermediate regimes. 

In all these cases, the light is produced in well-defined fiber modes, which is convenient for 
integration with waveguide devices. This holds the promise of combining many quantum logic 
gates using integrated optics, leading to the possibility of scalable quantum optical devices. 

2. Spontaneous four wave mixing theory 

In this paper we study SFWM in a single-mode optical fiber with a third-order nonlinear sus- 
ceptibility '. In this process, two-photons from pump fields E\ and Ei are jointly annihilated 
to create a photon pair comprised of one photon in the signal mode, E s , and one photon in 
the idler mode E(, We assume that all fields propagate in the fundamental spatial mode of the 
fiber. This assumption is justified if the fiber core radius is small enough that it only supports 
the fundamental mode, or alternatively if only this fundamental mode is excited. Following a 
standard perturbative approach |26|, the two-photon state produced by spontaneous four-wave 



mixing in an optical fiber of length L can be shown to be given by Il27ll28l 



|^> = |0),|0), + kJ J d(o s d(OiF (co s , a,) |a). s ),>,>, . (1) 

Here, K is a constant which represents the generation efficiency (linearly proportional to the 
fiber length, electric field amplitude for each of the pump fields and dependent on the relative 
polarizations of the pump and created pair fields) and F (o s , 0),) is the joint spectral amplitude 
function (JSA), which describes the spectral entanglement properties of the generated photon 
pair 
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which is given in terms of the pump spectral amplitudes OJi^ffl), and the phase mismatch 
function Ak (<Oi , <o s , (Oi), that in the case where the two pumps, signal and idler are co-polarized, 
is given by 

Ak (tai , C0 f , (Oi) =k((0i)+k((0 s + (0i-(0])-k((0 s )-k((0i)- (71 Pi + y 2 P 2 ) , (3) 

which includes a self/cross-phase modulation contribution for the two pumps with peak powers 
Pi and P2, characterized by the nonlinear parameters 71 and 72, which depend on specific fiber 
used and pump wavelength ||29l [30l . The energy conservation constraint is apparent in the 
argument of the second term of the phase mismatch (see Eq.([3]l). A "factorable" state is defined 
to be a state for which F (<O s , (Oi) is equal to a product of two functions, F (<O s , (Oi) = S((0 s )I((0i), 
where the functions S((o) and I((o) depend only on the signal and idler frequencies respectively. 

By making use of a linear approximation for the phase mismatch, in addition to modeling 
ai 2(ft)) as Gaussian functions with bandwidth Co respectively, it is possible to obtain an 
expression for the joint spectral amplitude in closed analytical form. Expanding fe(fi)^) in a 
first-order Taylor series about frequencies (0^ for which perfect phase-matching is attained 
(where ji — 1,2,5,/), and defining the detunings V s = <a s — (0® and v, = — (0®, the approximate 
phase mismatch Akn„ is defined by 



LAk lin = LAk^ +T s V s + TiVi, 



(4) 



where Afc(°), given by Eq.(|3j> evaluated at the frequencies 0)^, must vanish to guarantee phase 



matching at these center frequencies. The coefficients 7^ are given by = T ll + T p of/(o 
erf), where represent group-velocity mismatch terms between the pump centered at fre- 
quency 0), and the generated photon centered at the frequency 0)^, and z p is the group velocity 
mismatch between the two pumps 



T M =L 



z p =L 



^K)-4> 2 °) 



(5) 



written in terms of (<o) — d"kn /d<o" | m=co a . Note that this approach requires a priori knowl- 
edge, for given pump fields, of the signal and idler frequencies (0) 4 ° and G)f) at which perfect 





Fig. 1. (a) Pump envelope function a(co s ,(Qj) for a fiber characterized by r = 0.67/im, 
/ = 0.52 and L = 30cm. (b) Phasematching function (j) (v s , v,-); for a relatively small region 
of {co s , co,} space the phase-matching function contours are essentially straight lines, with 
slope s i = — arctan(r s /7)) (with 7^ given by Eq. (pi), (c) Resulting joint spectral intensity. 



phase-matching is achieved. These frequencies can be determined by solving (for example 
numerically) the condition Ak^ — (see Eq.^). It can be shown that within the linear ap- 
proximation, the integral in Eq. Q can be carried out analytically, yielding 



Fun (v s , V,) = a (V s , V;) (j) (V s , V/) . 



(6) 



where the pump envelope function a (v s , V,-) is derived from the pump spectral amplitudes for 
the two individual pump fields through the integral in Eq. and is given by 



a{v s ,Vi) = exp 



(v. s + v,f 



(7) 



and where (v s , V,) describes the phase-matching properties in the fiber. For degenerate pumps 
(where a\ (co) — 0.2(0))), it may be shown that the phase-matching function is given by 



<l>(v s ,Vi) = sinc 
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where LAku„ is given in Eq. (|4ji with T p = 0. For non-degenerate pumps, T p ^= and (y s , V;) 
4>(B;LAfc,,„), with 



erf ( —-iBx + erf (/fix) 
2B 



(9) 



where erf(z) is the error function, the parameter B is defined as B 



/(ai<T 2 T p ), 



and M is a normalization coefficient. In the present paper, we concentrate on the important class 
of factorable states for which F ((O s , 0),) = S (co s ) I (fflj) . 



3. Phase and group-velocity matching properties of photonic crystal fibers 



Our analysis focuses on PCFs consisting of a fused silica core surrounded by silica cladding 
with a pattern of air holes which remains constant along the fiber length. This mixture of air and 




Fig. 2. (a) Black, solid curve: phase-matching (Ak = 0) contour for SFWM in the degenerate 
pump case. Colored background: phase-matching orientation angle. Black, dashed line: 
symmetric group velocity matching (GVM) contour. Along the phase-matching contour 
we have indicated particular orientation angles of interest, (b) Close up, near the lower zero 
group velocity dispersion frequency. 



and glass in the cladding results in an average refractive index that is considerably lower than 
that of the core, providing a high dielectric contrast, resulting in strong optical confinement. 
This leads to high peak irradiances even for modest input powers, which enhances nonlinear 
optical effects such as SFWM. In addition, the dispersion characteristics of the PCFs can be 
engineered by variations of the distribution, size and shape of the air holes surrounding the 
core. In particular, it becomes possible to choose the zero dispersion wavelength(s) (ZDW), 
to tailor the FWM phase-matching properties [31 1 and to design fibers approaching endlessly 
single-mode behavior l32l . 

The FWM phase-matching properties are determined by the fundamental mode propagation 
constant, given in terms of the effective refractive index n e ff by k(co) — n e ff{(o)(o/c. We adopt 
a step-index model, where the core has radius r, its index is that of fused silica n s (co), and the 
cladding index is calculated as n c i al i(co) = / + (1 — f)n s (co), where / is the air-filling fraction. 
This fiber dispersion model has been shown to be accurate for / from 0.1 to 0.9 according to 
Ref. ll33l . In the context of our work, this model permits a straightforward exploration of the 
spectral entanglement properties in {r,f} parameter space. 

The structure of the SFWM two-photon state, in the degenerate-pump, co-polarized case (for 
which 71 = 72 = 7) is illustrated in Fig. [T] for a specific fiber with r = 0.67/xm, / = 0.52 and 
length L = 30cm. Fig. [TJa) shows the pump envelope function plotted as a function of <a s and 
COj (see Eq.(|7]i) where we have assumed that the pump is centered at 723nm, has a bandwidth 
of lnm, the incident power is 5W and 7= 70km _1 W _1 . Fig. |TJb) shows the phasematching 
function (see Eq.([8])) for this choice of parameters. Note that the phasematched region forms a 
strip on {co s , 0),} space, oriented at an angle VI - = — arctan(T s /T,) = —55° with respect to the 
(O s axis. Fig. |TJc) shows the resulting joint spectral intensity \Fn n ((D s ,(Oi)\ 2 (see Eq.jfrjl). It is 
apparent from Fig.[T]that the properties of the two photon state are determined by the i) relative 
orientations and ii) widths of the strips representing the phasematching and pump envelope 
functions. In this paper we explore how the interplay of the various design parameters can lead 
to two-photon states with engineered spectral entanglement properties. 



A plot of the perfect phase-matching contour (hk{co p ,(O s ,2(Op — (0 S ) = 0) for degenerate 
pumps versus co s and co p gives for each pump frequency the expected signal and idler central 
frequencies o 4 ? and G)f. Such a phase-matching contour is illustrated in Fig.|2]for r = 0.616/im 
and / = 0.6, with pump power 30W and y = 70km~ 1 W , where the generated frequencies 
are expressed as detunings from the pump frequency A s j — 0) 4 , — (O p . Note that energy con- 
servation implies that A s = — A,-. PCFs often have two ZDWs, where one can be as low as 500 
nm (in comparison with 1270 nm for bulk silica). The phase-matching contours take the form 
of closed loops, with inner branches near the pump frequency and outer branches that can be 
hundreds of nanometers from the pump wavelength. Note that four wave mixing relying on 
outer-branch phasematching has been observed in previous work, in the context of classical 
non-linear optics ifTSI . 

For this specific fiber geometry, a continuous pump wavelength range of approximately 453 
nm exists, in which parametric generation can be observed, almost completely contained be- 
tween the two ZDWs (X zl ii = 0.668/im and X z( /2 = 1.132/im). The power-induced phase modu- 
lation terms in Ak split the trivial A s = A, = branch, leading to frequency-distinct inner-branch 
signal/idler frequencies close to the pump (used for SFWM in iTPTl ). whereas the outer branch 
exhibits comparatively little dependence on pump power lfT5l[l6l[T7ll . 

In what follows we will concentrate on phase-matching configurations where the signal and 
idler frequencies are sufficiently removed from the pump frequency, or where the pump is or- 
thogonally polarized to the signal and idler photons, so as to avoid the generation of background 
photons by spontaneous Raman scattering at the lower-frequency photon of the pair [ 14 1. 

From Eq. (j7j) it is clear that the pump envelope function in {co s , ©,} space has contours of 
equal amplitude which have negative unit slope in all of {o. s , 0),} space. In contrast, the contours 
of the phase-matching function (co s ,COi) are characterized by a slope in {co s ,(Oi} space given 
by 9 S i — — arctan(r s /7}) (with given according to Eq. |5|l). The relationship between the 
slope 0Ap °f the curve at a point in {A s j,co p } space to the resulting phasematching function 
in {co s ,COj} space is 0& p = 45° — 9 S j. For example, 0„ = 45° corresponds to zero slope on 
the contour in Fig.|2ja) or 9a p = 0. The colored background in Fig. [^indicates the slope of the 
phase-matching contour, or orientation angle, ranging from 4 , = —90° to 0„ = +90°, indicated 
in blue and red respectively. 

The type of spectral correlations observed in a SFWM two-photon state is determined in 
part by the slope 0„ of the phase-matching contour. If the phase-matching contour is given 
by a closed loop (which is true for most fibers, i.e. range of values of {r,f}, of interest), all 
phase-matching orientation angles 9 s t are possible, controlled by the pump frequency. Thus, for 
certain relative orientations and widths of these two functions, it becomes possible to generate 
factorable two-photon states. It can be shown that a factorable state is possible if 



T s Tt < 0. 



(10) 



Among those states which fulfil Eq.dOl those exhibiting a phasematching angle of si = 
45°, or T s — —Tj, are of particular interest. For these states, in the degenerate pumps case, a 
factorable, symmetric state is guaranteed if 



2ra 2 \T s Ti\ = 1, (11) 

where T sa 0.193. The condition in Eq. ( 10 1 constrains the group velocities: either k^((O s ) < 
k^((Op) < ^((Oi), or k^((Oi) < k^((O p ) < k^((O s ) must be satisfied. The condition in Eq. 



(Hi constrains the bandwidths. Thus, the region in {o 4 , 0),} space in which factorability is 
possible is bounded by the conditions 7^ = and Tj = 0. 

Note that in the case of PDC (in second-order nonlinear crystals), it is in general terms chal- 
lenging to obey Eq. ( 10 1, for it implies that the pump, say with frequency co p , must propagate 



at a higher group velocity than one of the generated photons at (O p /2, which can be interpreted 
as anomalous group-velocity dispersion. In practice, this can be achieved for certain materials, 
within a restricted spectral region (usually with PDC in the infrared) in a type-II process, where 
the polarization of the pump is orthogonal to that of one of the generated photons |18j. In the 
case of SFWM, because the frequencies (in the degenerate pump regime) will in general obey 
(O s < (Op < (Oi or a>i < (Op < (0 S , no such anomalous group-velocity dispersion is needed, mak- 
ing it much more straightforward to fulfill the group-velocity matching conditions required for 
factorable photon pair generation. 

Once the dispersion relation for a given fiber geometry has been determined, the next step in 
designing a factorable photon pair source consists of identifying the pump wavelengths that sat- 
isfy the required group-velocity matching, amongst those which also satisfy phase-matching. 
For this purpose, we refer to the colored background in Fig. [2] Each color corresponds to a 
different phase-matching orientation angle 0„. In particular, symmetric factorable two-photon 
states, for which the signal and idler photons have identical spectral widths, are possible if 
8 S i = 45°. In this case symmetric group-velocity matching 2kp =ks + (or equivalently 
T s = —7}) is attained, and the phase-matching contours are oriented so that its contours have 
unit slope. The frequency values that fulfill this condition are represented by the dashed line. 
The pump frequency that permits symmetric factorable states can be determined from the in- 
tersection of the phase-matching contour with the group-velocity matching contour. 

Similarly, asymmetric factorable two-photon states, for which the signal and idler photons 
have greatly different spectral widths, are possible if J( = 0° or i7 = 90° . This is the case of 
asymmetric group-velocity matching, for which k p =ks or k p =kj (or equivalently T s = 
or Tj = 0) is required. In this case the phase-matching contours are oriented parallel to the (O s 
or (Oi axes. In addition, Eq. (Hi leads to the condition that 7} 3> l/<7 (for T s = 0) or T s 3> 1/(7 
(for Ti = 0). 

Both symmetric and asymmetric GVM are possible using a number of different configura- 
tions. We provide examples of each, exploiting co-polarized as well as cross-polarized SFWM, 
in the following sections. Note however, that the examples to be shown do not include all pos- 
sible source designs. 



4. Co-polarized fields and degenerate pumps: symmetric factorable states 

The synthesis of a symmetric, factorable state is illustrated in Fig.[3]for the same fiber geometry 
as described above (r = 0.616/xm and / = 0.6). Figure|3ja) shows the pump envelope function, 
for a single pump centered at 7 15 nm with a relatively narrow bandwidth of 0. lnm. Figure [3jb) 
shows the phase-matching function assuming a fiber length of 25cm, while Fig. j3jc) shows the 
joint spectral intensity \F[j n ((O s ,(Oj)\ 2 , exhibiting an essentially factorable character. For com- 
parison, Fig. [3jd) shows the joint spectral intensity obtained by numerical integration of Eq. 
Q using the full dispersion (rather than relying on a linear approximation), revealing that this 
case the linear approximation of the phase-matching condition is in fact an excellent approx- 
imation. This source leads to a numerically-obtained state purity defined as Tr[p : 2 ], where p v 
is the reduced density operator for the signal state, of 0.901. In fact, the departure from ideal 
purity is mainly due to sidelobes (related to the sine function) displaced from the central portion 
of joint spectral intensity shown in Figs. [3|c) and (d). The purity can be increased by filtering 
out these sidelobes, which in general contain a small fraction of the total flux. Specifically, for 
the example shown in Fig. [3] two separate narrowband rectangular-profile spectral filters for 
the signal and idler modes with equal frequency bandwidth of 9.35 x 10 n rad s (correspond- 
ing to wavelength widths AA S « O.lnm and AA, w 1.2nm) increases the purity to 0.981, while 
reducing the flux by < 5%. 

As pointed out already, one of the advantages of generating factorable, symmetric photon 
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Fig. 3. Joint spectral intensity (JSI) obtained for the fiber geometry assumed for Fig. [2] 
(r = 0.616/lm and / = 0.6) where the pump central wavelength (X p o = 0.7147/im) is ob- 
tained by imposing simultaneous phase-matching and group-velocity matching. We con- 
sider a pump bandwidth of 0.1 nm and a fiber length of 0.25m. The function values are 
normalized such that white = 1 and black = 0. (a) Pump envelope function cc(co s , CO,-), (b) 
Phase-matching function <j> (v s , V/). (c) Analytic JSI, obtained with approximation from Eq. 
p). (d) JSI obtained by numerical integration of Eq. |2| . 

pairs through SFWM rather than PDC in a second-order nonlinear crystal, is that it presents 
greater (but still not arbitrary) choice of the pump, signal and idler frequencies. Thus, for the 
example shown above, we have chosen the pump so that it is compatible with, say, a Ti:sapphire 
laser. One of the generated photons is centered at 1551 nm, which corresponds to the stan- 
dard telecommunications band and could be detected with an InGaAs-based avalanche photo- 
detector. The conjugate generated photon is centered at 464 nm, which can be detected with a 
silicon-based avalanche photodiode (though with non-optimal efficiency, since this wavelength 
is relatively far from the sensitivity peak of such detectors). A natural application for such a 
source would be the heralding of pure, single photons for transmission over standard optical 
telecommunication fibers. 

5. Co-polarized fields and non-degenerate pumps: symmetric factorable states 

The use of non-degenerate pumps represents an additional degree of freedom, not available for 
PDC in second-order crystals, which can be exploited for state engineering purposes. The two 
pumps may differ both in terms of their central frequencies (a® and CO®) and their bandwidths 



(ffi and (J2); in what follows we exploit both of these aspects of non-degeneracy. The orientation 
of the phase-matching function (see Eqs. Q and (|9|) in {© v ,0),} space is determined by the 
angle S , = — arctan (T s /Tj). Note that the orientation of the phase-matching function depends 
on oj and 02 (through T s and 7)), while the shape of its profile and its width also depend 
on (7i and 02 (through parameter B). A considerable simplification results by imposing the 
condition G\ <C 02- In this case, 7Jj reduces to the corresponding degenerate-pump values 
(with jj. = s,i). Furthermore, parameter B reduces to l/(oiT p ), so that the phase-matching 
function width exhibits no dependence on the broader pump bandwidth (02), while the effective 
pump envelope function depends only on 02- This de-coupling of (Jj and 02 translates into a 
more straightforward exploration of possible fiber geometries for the generation of states with 
specific spectral entanglement properties. 

A phase-matching diagram similar to that presented for degenerate pumps (see Pig. |2j, may 
be prepared by maintaining one of the two pump frequencies fixed (in this case pump 2), where 
it is convenient to express the generated frequencies as detunings from the mean pump fre- 
quency i.e. A s ; = (0 s .i — (oj 1 + Gk)/2; note that energy conservation implies that A s = —A;. 
Figure|4]illustrates the phase-matching properties for a specific non-degenerate geometry, with 
r = 0.601^m, / = 0.522 and a% = 1.508 x 10 15 rad/sec (that is, A 2 ° = 2nc/o^ = 1250 nm). 
The black solid curve represents the resulting phase-matching contour. The two straight lines 
represent trivial phase-matching branches (in which the created photons are degenerate with 
the two pumps), while the loop represents the non-trivial phase-matching branch. Depending 
on the fiber geometry, the non-trivial branch may become large enough to overlap the trivial 
branches, or may shrink down to a single point. For the former case, if self-phase modulation 
becomes appreciable, this contour splits into three distinct loops (rather than two as for degen- 
erate pumps). The colored background represents the orientation angle 8 S i ranging from —90° 
in blue to 90° in red. 

A factorable, symmetric two-photon state is produced if i) 9 S i = 45° which implies T s = —Tj 
and ii) the spectral widths of the phase-matching function (j) (co s , o,-) and the pump envelope 
function a((D S: fflj) are equal to each other. From Eqs. ^ it is straightforward to show that the 
first condition is satisfied if the four fields satisfy 

Ik^-k^-k^ = -2{k^-k^)a 2 J{al + al). (12) 

This represents a generalization of the equivalent group-velocity matching condition for de- 
generate pumps (see Eq.flO] and the text which follows). If <Ti <C 02, the right-hand side of Eq. 



(12 1 vanishes, leading to a condition that is identical in form to that obtained for degenerate 
pumps (where the broader pump now plays the role of the degenerate pump). In Fig. |4] pairs 
{A,, ©pi} that satisfy this generalized group-velocity matching condition are represented by a 
dashed line. Note that the position and shape of this contour in general depend on <3\ and (72; 
however, in the case <3\ -c 02, the contour becomes decoupled from these bandwidths. Inter- 
section points between the two contours determine the center frequency for pump field 1 that 
satisfies the requisite group-velocity matching, in addition to phase-matching. The values of r 
and / are selected so that 0„ = 45° is satisfied and so that the broadband pump frequency (0® 
takes a certain desired value (in this case, — 2nc/(0^ — 0.625jtim). This leads to the gen- 
erated signal and idler wavelengths 736nm and 960 nm. The fiber length, along with the two 
bandwidths (subject to Ci <C 02) are selected so as to match the widths of the pump envelope 
and phasematching functions. For the specific geometry shown, these values are L = 25cm, 
while (Ji and 02 are given in terms of the corresponding FWHM bandwidths: AAi =1.51 nm 
and AA2 = 0.12 nm. 

An important motivation for exploiting pump non-degeneracy is that it permits SFWM ge- 
ometries with pumps which are sufficiently distinct from the signal and idler photons, and yet 
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Fig. 4. Black, solid curve: phase-matching (Ak = 0) contour for SFWM in the non- 
degenerate pump regime. Colored background: phase-matching orientation angle. Black, 
dotted line: frequencies that satisfy the group-velocity matching condition (see Eq. j!2[>). 
Along the phase-matching contour we have indicated particular angles of orientation of 
interest. 

where the latter are in relative proximity. This is important on the one hand for the suppression 
of contamination from spontaneous Raman scattering and on the other hand to ensure that both 
generated photons can be detected with available efficient single-photon detectors. The pump 
wavelengths were selected so that they can be obtained from a modelocked Cr:Forsterite laser 
and its second harmonic, while the generated light is within the detection bandwidth of sili- 
con avalanche photo-diodes. Figure |5|a) shows the pump envelope function <x((O s , (fy). Figure 
|5|b) shows the phase-matching function (co s , Cfy) (see Eq.(|9|) for which B = 1 .73. Figure^c) 
shows the resulting joint spectral intensity which exhibits a factorable character, obtained using 
the linear-dispersion approximation (see Eq. Q). For comparison, Fig. [5jd) shows the joint 
spectral intensity obtained by numerical integration of Eq.Q, exhibiting excellent agreement 
with Fig. |5jc). This source leads to a numerically-obtained state purity of 0.89. As in the case 
of a factorable symmetric state obtained with degenerate pumps (see Fig. [3}, the purity can be 
increased by filtering out the sidelobes in the joint spectral intensity at small cost in terms of 
collected flux. Two separate narrowband rectangular-profile spectral filters for the signal and 
idler modes with equal frequency bandwidth of 12.90 x 10 12 rad s^ 1 (corresponding to wave- 
length widths AAy w 3.7nm and AA, « 6.3nm) increases the purity to 0.98, while reducing the 
flux by < 5.9%. 
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Fig. 5. Joint spectral intensity (JSI) \F(co s , ft>;)| 2 obtained for the fiber geometry assumed 
for Fig. [5] (a) Pump envelope function a(co s , caf). (b) Phase-matching function 0(tO s ,<Bj). 
(c) Analytic JSI. (d) JSI obtained by numerical integration of Eq. |2j. 

6. Cross-polarized fields and degenerate pumps: asymmetric factorable states 

Birefringence is an additional degree of freedom that we can use to satisfy the conditions nec- 
essary for generating factorable two-photon states. Typical hexagonal pattern PCF can exhibit 
significant birefringence due to slight asymmetries in the hole pattern or internal stresses in the 
glass. Birefringence can also be purposefully added to the structure through defects such as two 
larger holes or elliptical cores. In PCFs, the birefringence typically ranges from An w 1 x 10~ 5 , 
for normal PCF, to An f=s 7 x 10~ 3 |34| for polarization maintaining PCFs. 

Conservation of angular momentum restricts the types of FWM processes that can occur in a 
fiber, with the result that the polarizations of the four fields must appear in pairs in the process 
||29ll35ll . It follows that there are three distinct pair generation processes that can occur in fibers: 
xx — * xx, xy — > xy, xx — > yy, where x and y can be any two linear orthogonal polarizations. Photon 
pair generation with frequency separated pumps in the xy — > xy process was analyzed in 2004 
11361 and demonstrated [9 1 in 2005. In contrast, we limit our discussion to a single pump pulse 
polarized along one fiber axis (e.g. xx — > yy, 2(O p — > 0), + CO s ). The first realization of birefringent 
FWM was a demonstration of this process in regular fibers |37l . It was demonstrated in PCFs 
in 2006(381 • 

Single-photon detectors operating in the visible region of the spectrum are generally less 
costly and have higher quantum efficiencies and lower noise than single-photon detectors oper- 



ating in the near-infrared region. Therefore for some applications it is desirable for both photons 
from a SFWM pair to be produced in the visible region. This limits the usefulness for these ap- 
plications of PCFs which have outer-branch solutions far separated from the pump frequency, 
as shown in Fig. |2ja), as one of the photons is generated in the infrared. On the other hand, the 
photons should be at least 50 THz from the pump to reduce background from spontaneous Ra- 
man scattering. This limits the usefulness of the inner branch solutions, for which the generated 
photons are typically less then 5 THZ from the pump 1131 . The benefit of birefringence is that 
it creates new factorable solutions inside the visible region yet still separated from the pump 
by more than the Raman bandwidth. Unlike the nondegenerate case, it accomplishes this with- 
out the experimental complexity of two pumps. This can be seen by considering the relevant 
phase-matching equation, 



2 

Ak (03 P ,A S ) = 2k x ((o p ) - k y ((o p + A,.) - ky ((Op - A s ) --yP (13) 
m 2k x ((O p ) - k x ((O p + A,.) - k x ((O p - A s ) +2An°^- ^yP, (14) 



where the power-induced term is 1/3 of that in the co-polarized case [29 1 . Some insight can 
be gained by separating the birefringent contribution to Ak from the individual wavevectors as 



in Eq. (14i. It was experimentally observed in [38| that An is approximately independent of 
frequency and we note that (O p is relatively constant in comparison to k x ((O p ). Consequently, 
the birefringent contribution to the wavevector mismatch functions much like the power phase- 
modulation term in Eq. ([3]), introducing in the Ak = curve a splitting of the A s = solution 
(creating the inner branch in Fig . [2J . Unlike the power term, the birefringence can create a large 
splitting and can be either positive or negative. 

We begin by extending the step index model reviewed in Section 3 to birefringent fiber. In 
this extension, each transverse axis of the fiber is modeled as a separate fiber. The axes have the 
same air-filling fraction parameter / but different core diameters d and d + Ad, simulating the 
birefringence [33 , 38 1. We use this model for all the figures in this section. We also neglect the 
power term in Eq. ( p"3j ) since it is typically much smaller than the birefringent contribution. Al- 
though this model is accurate only for small birefringence where the fiber asymmetry is small, 
it allows us to predict universal properties of the SFWM across a broad range of configurations. 

In Fig. [6] the phase-matching (Ak((0p,A s j) — 0) and group-velocity matching (T s , = 0) 
curves are plotted for a PCF fiber with ZDWs at 790 nm and 1404 nm, where d = 1.75/xm 
and / = 0.43. Three pairs of curves are plotted for Ad = —0.001,0,0.001 fim, resulting in a 
birefringence of An = -3 x 10~ 5 ,0,3 x 10~ 5 (where An = riy — n x ). The corresponding shift in 

the ZDW is -0.07, 0, 0.07 nm and k^ —ky — —0.5,0,0.5 ps/km, respectively. The segments of 
the Ak = curve bounded by intersections with the T s j = curve are regions where factorability 
is possible. For both positive and negative An the outer branch intersections remain approxi- 
mately the same as in the copolarized degenerate case. In addition, negative An creates two 
intersections (Points C and D) in between the two ZDWs. In between these points factorability 
is possible. Likewise, positive An also creates two intersections (Points E and F). However, now 
it is outside the region bounded by the intersections that factorability is possible. The splitting 
of the phase-matching (Ak = 0) solution from the A s j — line in the An < case is similar to 
that from the power-induced phase modulation term 2yP in the phase matching equation Eq. 
Q, but is an order of magnitude larger. Unlike power-induced inner branch FWM, birefrin- 
gence allows for phase-matching beyond the Raman peak and is relatively insensitive to pump 
laser power fluctuations. In contrast, the An > case is qualitatively different from either of 
these since 2yP is always positive in silica. Moreover, it allows for factorable state generation 
over an unprecedented pump wavelength range, essentially anywhere outside the ZDWs. 




Fig. 6. a) Cross-polarization phase-matching Ak = (thick) and group-velocity matching 
T sJ = (thin) curves ( / = 0.43, d = 1 .75/Xm, ZDWs = 790 nm, 1404 nm - similar to NL- 
1.7-790 from Crystal Fiber). With power P = 0, three pairs of curves are plotted with Ad = 
—0.001,0,0.001 yon, respectively resulting in a birefringence of An = — 3 x 10~ 5 ,0,3 x 
10~ 5 (red, green, blue). Points A through F bound regions in which factorability is possible, 
b) Full dispersion numerical JSIs for the asymmetric states corresponding to points D (0 S! - = 
0° , A s = 746.4 nm and A; = 949.0 nm, AX P = 0.30 nm, and L = 30 m) and F (G si = 90°, A. s 
= 677.1 nm, A; = 834.2 nm, AA P = 0.30 nm, and L = 30 m). 



Figure |6jb) shows the joint spectral intensity in two cases that lead to high degrees of fac- 
torability. Factorability corresponds in both cases to the absence of correlation between signal 
and idler frequencies, as can be seen in the shape of the JSI. In either case, if the idler photon 
is detected (without spectral filtering), then a nearly pure-state signal photon is heralded. Nu- 
merically calculated heralded photon purity for the cases at point F and point D are 0.988 and 
0.993, respectively. 

In Fig. |7J the intersection of the Ak = and T s j = curves is plotted as a function of the 
An at 800 nm. Thus for each An the figure indicates the pump wavelength X p that creates an 
asymmetric factorable state (9 s i = 0° or 8 s i = 90°), as well as the center wavelengths A s and A; of 
the generated photons. At these points, the pump bandwidth must satisfy condition T v , >> a -1 
but is otherwise arbitrary. The shaded regions indicate the X p range in which factorability is 
possible (including the symmetric state) if the pump bandwidth is set according to Eq. (Hi. 

It has been pointed out in the FWM literature that birefringent FWM has the added benefit 
of reducing Raman contamination in the signal and idler [39 1. Recently this has also been 
forwarded as a strategy to reduce background for pair-generation BOl . Raman gain gRaman is 
reduced by as much as an order of magnitude in the axis orthogonal to the strong pump pulse 
lfl4l . However, the relevant %^ for birefringent pair production is also reduced by a factor 
of three compared to co-polarized FWM [29|. Tripling the power P to compensate, Ref. ATI 
studied the pair-correlation produced by the fiber, p c (0) = jp^pj — 1 , where I is the intensity of 
signal or idler mode. They found that at the Raman peak, p £ (0) was 7 for co-polarized SFWM 
as compared to 60 for cross-polarized SFWM, indicating that the cross-polarized case can be a 
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Fig. 7. Points C through F from Fig. [6] are plotted as a function of the birefringence An. 
Thus, the solid dark blue curves give the pump wavelength X p of asymmetric factorable so- 
lutions. The corresponding idler and signal wavelengths A; and X s are given by the magenta 
dashed curves (for points D and F) and light-blue dash-dotted curves (for points C and E). 
The shaded regions indicate the range in which factorable states can be generated provide 
the pump bandwidth is matched to the fiber length. The fiber parameters are the same as in 
Fig.© 

higher quality photon-pair source. 

7. Co-polarized fields: ultra-broadband two-photon states 

The flexibility afforded by PCFs in engineering the spectral correlations of photon pairs permits 
the generation of a wide class of states. For instance, this may be illustrated by considering the 
generation of quantum states in which the photon pair component has an extraordinarily high 
degree of frequency anti-correlation[42|. The photons individually have an ultrabroad band- 
width, yet are very tightly anticorrelated in frequency, and highly correlated in time (to the level 
of a single optical cycle). We have identified conditions for which ultra-broadband two-photon 
states can be generated in the degenerate pump regime by SFWM. Such ultrabroad states have 
recently been generated by PDC ll43l . As pointed out above, for fiber geometries where the 
phase-matching diagram occurs in the form of loops (which is true for most cases of interest), 
it is possible to obtain an arbitrary phase-matching orientation on {co s , ©,} space. By imposing 
constraints on higher dispersive orders, it becomes possible to generate states with yet different 
character. In particular, in analogy with the generation of ultra-broadband two-photon states 
via PDC l43ll44l . it can be shown that in the case of SFWM, the bandwidth of the generated 
light can be made particularly large if the fourth-order and second order coefficients 
vanish simultaneously. In this case, a degenerate pump at the zero-group-velocity-dispersion 
frequency can lead to a state with a remarkably broad bandwidth (hundreds of nanometers), 
even for a nearly monochromatic pump. For example, for the PCF described previously with 
r = 0.616/im and / = 0.6, whose phase-matching curve is shown in Fig. [2] pumping at a fre- 
quency around 2.82 x 10 15 rad/sec (wavelength A = 668 nm), where the curve is nearly vertical, 



would create photons with a bandwidth of around 2 x 10 15 rad/sec (AA » 542 nm) centered at 
the pump frequency, depending on the bandwidth of the pump. Previously, light produced by 
co-polarized pumps was shown to have a bandwidth of 50 nm for CW pumps [45 1 and 100 nm 
for pulsed pumps |46|. For cross-polarized CW pumps broadband light of 34 nm was produced 
ED- 

8. Conclusions 

We have studied theoretically the spontaneous four wave mixing process in single mode pho- 
tonic crystal fiber. We have shown that it is possible to design a source so as to yield photon 
pairs with a broad class of spectrally engineered properties, including factorable and ultra- 
broadband states. These results are achieved by adjusting the source configuration, including 
the dispersion (determined by the core radius and air filling fraction) and birefringence of the 
fiber and the pump frequency or frequencies, bandwidth(s) and polarization(s). The flexibility 
of the two-photon state engineering arises from the interplay of this large number of control pa- 
rameters. Factorable signal-idler pair generation allows the production of heralded, high purity 
single-photon wavepackets, with minimal or no need for spectral filtering. 

Factorable state generation is possible if certain group-velocity matching conditions are 
satisfied, in addition to standard phase-matching. Most fibers (all those for which the phase- 
matching curve occurs in loops) can achieve all appropriate group- velocity matching conditions 
by suitable choice of the pump wavelength. This allows one to realize all possible orientations 
of the joint spectral amplitude of the generated photon pairs. The size of the phase-matching 
loop is largely determined by the dispersion of the fiber, hence the utility of PCFs. Therefore 
we find that for specific choices of pump frequency most fiber designs can yield factorable 
states. This is in contrast with parametric downconversion in second-order nonlinear crystals 
where such conditions are met at best at isolated frequencies for certain materials, or using 
more complicated microstructures Il4l l23~l . 

An important design consideration is the possible contamination of the frequency sideband 
modes by photons from spontaneous Raman scattering. This problem can be circumvented in 
part by cooling the fiber to cryogenic temperatures, thereby suppressing the Raman scatte- 
ring iflOl [481 . We have shown how to circumvent the problem by exploiting a phase-matching 
configuration where the signal and the idler frequencies are widely separated from the pump 
frequency (see Fig. |2j. If desired, the phase matching in non-birefringent fiber can be engi- 
neered to create photon pairs having one photon in the infrared and one photon in the visible. 
Alternative approaches also discussed in this paper are cross-polarized FWM in birefringent 
fiber, where the pump and generated photons have orthogonal polarizations, and the use of 
non-degenerate pumps in a co-polarized regime. We have shown that these approaches lead to 
sources capable of generating factorable photon pairs where the signal and idler photons are 
both in the visible, and at the same time are both sufficiently separated from the pump fre- 
quency to avoid Raman contamination. We expect these results to be useful in the design of 
sources for practical implementations of quantum information processing technologies. 
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